DDC  FILE  COPY  AD  AO  58818 


\q/— — \>  . „ 

[AFOSRjfRr^l-^-^  3 6 Z- 


AFOSR  Scientific  Report 


MAE-TR-198 


LEVEL 


-^%AN  ANALYTICAL J10DEL  FOR  BONDED 

/joint  analysis  in  Composite  Structures 

INCLUDING  HYGROTHERMAL  EFFECTS. 


Q/J ZtitzviM  re?*'; 

(fo)^ Robert  C. [Wetherhold 

/ ,Jack  R.  Vinson 


Department  of  Mechanical  and  Aerospace  Engineering 
University  of  Delaware 
Newark,  Delaware  19711 


\U  Febiwuwy  H77 


f; 1 


Air  Force  Office  of  Scientific  Research 

United  States  Air _Eqrc£ -7 

Contract  N>h^h>/AF0SR^7A-2759  / 

Approved  for  public  release;  distribution  unlimited 


43 


* 


4b'  ■- 

Qualified  requestors  may  obtain  additional  copies 
from  the  Defense  Documentation  Center,  all  others 
should  apply  to  the  National  Technical  Information 
Service. 


Conditions  of  Reproduction 


Reproduction,  translation,  publication,  use  and  disposal 
in  whole  or  in  part  by  or  for  the  United  States  Govern- 
ment is  permitted. 


AIR  FORCE  OFFICE  OP  SCIENTIFIC  RESEARCH  (AFSC) 
NOTICE  OF  TRANSMITTAL  TO  DDC 
This  teohnioal  report  has  besu  reviewed  and  Is 
approved  for  put  1 <-«.•  r-'Iease  IA<V  APR  190-12  (7b). 
Distribution  is  unlimited. 

A.  D.  BLOSE 

Technical  Information  Officer 


^classified 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (When  Dele  l 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

i report  number 

AFOSR-TR-  7 8 - 13  8 7 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

AN  ANALYTICAL  MODEL  FOR  BONDED  JOINT  ANALYSIS  IN 
COMPOSITE  STRUCTURES  INCLUDING  HYGROTHERMAL 

EFFECTS 

5 TYPE  OF  REPORT  & PERIOD  COVERED 

INTERIM 

6 PERFORMING  ORG.  REPORT  NUMBER 

7 AuThOROJ 

ROBERT  C WETHERHOLD 

JACK  R VINSON 

8.  CONTRACT  OR  GRANT  NUMBERCfJ 

AFOSR  74-2739  ^ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

UNIVERSITY  OF  DELAWARE 

MECHANICAL  & AEROSPACE  ENGINEERING  DEPARTMENT  " 
NEWARK,  DELAWARE  19711 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 

2307B1 

61102F 

It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH/NA 

BLDG  410 

BOLLING  AIR  FORCE  BASE,  D C 20332 

12.  REPORT  DATE 

Feb  77 

13.  NUMBER  OF  PAGES 

13$ 

'*  MONITORING  AGENCY  NAME  A AOORESSf/f  dlllerent  Irom  Controlling  Otlice ) 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

15a.  DECLASSIFICATION  DOWNGRADING 

schedule 

16  DISTRIBUTION  STATEMENT  (oi  this  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  If  different  from  Report) 

i8  supplementary  notes 

19  KEY  WORDS  (Continue  on  reverse  side  it  neceaaery  end  identity  by  block  number ) 

COMPOSITE  MATERIALS 
BONDED  ADHESIVE  JOINTS 
LAMINATED  STRUCTURES 


0 ABSTRACT  f Continue  on  reverse  aide  if  necessary  and  identify  by  block  number) 

A plane-strain  model  is  developed  for  symmetric,  anisotropic  laminated  plates 
subjected  to  thermal,  moisture,  and  mechanical  loading  including  surface  tractionjs 
Provision  is  made  to  allow  assessment  of  the  relative  importance  of  the  normal 
strain  and  of  various  off-axis  effects  by  the  use  of  tracing  constants.  The 
solution  for  stresses  and  displacements  is  given  for  a single  plate  element,  and 
the  computations  are  demonstrated  by  numerical  examples.  The  basic  equations  for] 
a bonded  joint  composed  of  pl.me-st  ra  in  plate  adherends  are  derived  and  discussed.[ 
References  to  alternate  derivations  are  given  whenever  appropriate  in  order  to 
pp  form  H73  I facilitate  comparisons.  | 


SECU 


ASSl 


UNCLASS  I 

FlCATlON  OF 


SJHED 


THIS  PAGE  (When  Data  Entered) 


AN  ANALYTICAL  MODEL  FOR  BONDED 
JOINT  ANALYSIS  IN  COMPOSITE  STRUCTURES 
INCLUDING  HYGROTHERMAL  EFFECTS* 


CENTER  FOR  COMPOSITE  MATERIALS 
AND  THE 

DEPARTMENT  OF  MECHANICAL  AND  AEROSPACE  ENGINEERING  ^ 
UNIVERSITY  OF  DELAWARE 


February  1977 


* Research  sponsored  by  the  Air  Force  Office  of  Scientific 
Research,  Office  of  Aerospace  Research,  United  States 
Air  Force,  Grant  Number  AFOSR-74-2739  / This  document 
has  been  approved  for  public  release;  distribution  is 
unlimited. 

**  Research  Assistant 


***  Professor  and  Chairman 


I 


TABLE  OF  CONTENTS 


Abstract 

Nomenclature 

Figures 

CHAPTER  1 Introduction  and  Summary  of  Basic 
Assumptions 

CHAPTER  2 Derivation  of  the  Governing  Equations 
for  a Plane-Strain  Plate  Element 

2.1  Basic  Elastic  Relations  and  Discussion  of 
Plane-Strain  Condition 

2.2  Equilibrium  Equations 

2.3  Presumed  Form  for  Transverse  Shear  Stress 
Derivation  of  Normal  Stress  o 

z 

2.4  Calculation  of  In-Plane  Displacement,  u, 
and  Local  Slope,  3w/3x 

2.5  Derivation  of  the  Integrated  Axial  Load, 
Shear,  and  Moment  Resultants 

2.6  Derivation  of  Equilibrium  Equations  for 
Integrated  Resultants 

2.7  Coefficients  for  a Specific  Hygrothermal 
Case 

2.8  Summary  of  Formulae  for  Stress;  Calculation 
of  Strains 


CHAPTER  3 Studies  of  a Single  Plate  Element 

3.1  Solution  for  a Single  Plate  Element  Subject 
to  Known  Surface  Stresses 

3.2  Parametric  Studies  of  a Single  Plate  Element 
Subject  to  Known  Surface  Tractions  and  Known 
Temperature  and  Moisture  Profiles 


IV 


xz 


Page 

vi 

vii 

x 

1 

4 

4 

13 

14 

21 

52 

68 

71 

77 

79 

79 

87 


j 


CHAPTER  4 Derivation  of  the  Governing  Equa- 
tions for  a Bonded  Joint 


4.1  Development  of  Alternate  Expressions 
Employing  au 

4.2  Derivation  of  Interface  Equations  in 
Bonded  Joints 

4.3  Discussion  of  Bonded  Joint  Equations 
CHAPTER  5 Conclusions 

REFERENCES 

APPENDIX 

Comparison  of  this  research  with  the 
previous  analysis  by  W.J.  Renton  and 
J.R.  Vinson  as  presented  in  AFOSR  TR 
75  0125,  Aug.  1974 


ACCESSION  for  | 

MIIS  White  Section  ^ 

DDC  Buff  Section  □ 

t'NAWJNCrn 
JUST  *CAi  ion 


RY 

GIOrattMJIT  G02ES 


ABSTRACT 


A plane-strain  model  is  developed  for  symmetric, 
anisotropic  laminated  plates  subjected  to  thermal, 
moisture,  and  mechanical  loadings  including  surface  trac- 
tions. Provision  is  made  to  allow  assessment  of  the  rela- 
tive importance  of  the  normal  strain  and  of  various  off- 
axis  effects  by  the  use  of  tracing  constants.  The  solu- 
tion for  stresses  and  displacements  is  given  for  a single 
plate  element,  and  the  computations  are  demonstrated  by 
numerical  examples.  The  basic  equations  for  a bonded 
joint  composed  of  plane -strain  plate  adherends  are  derived 
and  discussed.  References  to  alternate  derivations  are 
given  whenever  appropriate  in  order  to  facilitate  compari- 
sons . 
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NOMENCLATURE 


k subscript 


m superscript 


total  thickness  of  the  laminate, 
z distance  from  the  geometrical  mid- 
plane of  the  laminate  to  the  inter- 
face between  lamina  (K)  and  lamina 
(K+l) . 

th 

denotes  the  k lamina, 
percent  of  moisture  content  measured 
from  some  stress  free  condition;  a 
smooth  function  of  x and  z. 
denotes  moisture  dependent  term, 
terms  of  the  lamina  stiffness  matrix 
referred  to  the  principal  material  axes, 
terms  of  the  lamina  stiffness  matrix 
referred  to  an  arbitrary  set  of  axes, 
terms  of  the  lamina  compliance  matrix 
referred  to  the  principal  material 


coordinate  system. 

terms  of  the  lamina  compliance  matrix 


referred  to  an  arbitrary  set  of  axes, 
transformation  matrix  used  between  the 


material  principal  coordinates  and  the 
geometric  coordinates. 


Vlll 


T Superscript 
u° 


w° 


1'  2'  3 


a , a , a , 
x y z 


xy 


K 


<V  S3 


ex.  ey.  bz, 

8xy 


AT 


denotes  temperature  dependent  term, 
the  laminate  mid-plane  longitudinal 
displacement  in  the  x direction, 
the  laminate  mid-plane  transverse 
displacement  in  the  z direction, 
coordinate  in  the  thickness  direction, 
measured  from  the  geometric  midsurface, 
coefficients  of  thermal  expansion, 
along  principal  material  axes, 
coefficients  of  thermal  expansion 
with  respect  to  laminate  coordinate 
system. 

boundary  constant  for  continuity  of 

o between  laminae, 
z 

coefficients  of  moisture  or  hygroscopic 
expansion,  in  material  principal  axes, 
coefficients  of  moisture  or  hygroscopic 
expansion  with  respect  to  laminate 
coordinate  system. 

temperature  measured  from  some  stress 
free  condition;  a smooth  function  of 
x and  z . 

normal  dilitational  strains 
shear  strains,  (i^j) . 
normal  stresses, 
shear  stresses,  (i^j). 


FIGURE  I 


Specimen  and  Fiber  Coordinate  System 
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T. 

T 1 | t 


FIGURE  2 

Integrated  Stress  Resultant  and 
Surface  Traction  Definitions 


ADHESIVE 


FIGURE  4 


Lamina  and  Height  Definitions  for  a Symmetric 
Plate  with  an  Odd  Number  of  Laminae 
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Combine  these  layers,  then 
proceed  as  in  figure  4 


FIGURE  5 

Symmetric  Plate,  Even  Number  of  Laminae 
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CHAPTER  1 


Introduction  and  Summary  of  Ba,sic  Assumptions 

In  many  advanced  structural  applications,  the  ad- 
vantages of  high  strength-to  weight  fibrous  composites  are 
well-proven.  What  is  often  lacking,  however,  is  a suitable 
method  of  analysis  for  laminated  composite  plates  bonded 
into  a structure  and  subjected  to  a full  range  of  mechanical, 
thermal,  and  moisture  loadings.  The  usual  plane  stress  plate 

theory  (t  = a_  = x _ = 0)  is  not  suitable  because  of  its 
x z z yz 

failure  to  account  for  the  surface  tractions  produced  by 
bonding.  The  theory  developed  herein  is  based  on  a plane- 
strain  approach,  which  would  be  suitable  for  many  panel 
constructions.  [Ref.  1,2,3] 

If  the  composite  is  polymeric,  thermal  and  mois- 
ture effects  (termed  "hygroscopic")  become  important. 

The  magnitude  of  the  induced  strains  due  to  these  effects 
as  well  as  material  property  degradation  at  elevated  tem- 
perature demand  their  inclusion  in  the  analysis.  A general 
procedure  has  thus  been  developed  to  predict  stresses, 
displacements,  and  strains  for  the  symmetric  laminated  plate 
under  combined  thermal,  moisture,  and  mechanical  loadings. 


2 


In  order  to  obtain  a solution  to  the  laminated 
plate  problem,  certain  assumptions  have  been  found  useful. 
These  assumptions  are  summarized  by  the  following: 

1)  A plane-strain  condition  is  assumed  in  the  y direc- 
tion, and  y dependency  of  all  displacements  and  of  all 
stresses  is  assumed  zero.  (See  Figure  1)  The  three  di- 
mensional elasticity  problem  is  thus  reduced  to  two  dimen- 
sions . 

2)  The  laminate  is  assumed  symmetric  with  respect  to  its 
midplane;  however,  either  an  odd  or  an  even  number  of  layers 
may  be  accomodated.  (See  Figures  4 and  5) . 

3)  The  material  coordinate  system  (1,2,3)  is  a Cartesian 
system  with  3 mutually  orthogonal  planes  of  symmetry. 

4)  The  transverse  shear  stress  distribution  in  the  lami- 
nate may  be  represented  by  a quadratic  function  in  z 
times  a function  of  x . (See  Figure  6) . 

5)  The  theory  presented  here  is  linear;  the  results  of 
in-plane  loadings,  surface  tractions,  and  thermal  and  mois- 
ture loadings  may  therefore  be  superimposed. 

6)  The  adhesive  between  the  plate  elements  is  subject  to 
shear  and  normal  stresses  only  as  shown  in  Figures  3 and  16. 


i 


3 

7)  Both  transverse  shear  strain,  e , and  normal  strain, 

x z 

ez,  are  included.  Provision  is  made  to  neglect  the  normal 
strain,  if  desired,  by  use  of  a tracing  constant. 

8)  The  given  temperature  and  moisture  functions,  AT 

and  m,  will  be  assumed  smooth  (possessing  continuous  deri- 
vatives) in  x,  y and  z throughout  the  laminate,  and 
also  will  be  of  polynomial  form  in  z . 

9)  All  material  properties  are  assumed  constant.  Should 

! 


high  temperature  or  moisture  alter  the  properties  signifi- 
cantly, average  properties  will  be  used. 


CHAPTER  2 


Derivation  of  the  Governing  Equations  for 
a Plane-Strain  Plate  Element 


2.1  Basic  Elastic  Relations  and  Discussion  of  the  Plane- 
Strain  Condition 


Given  a constant  coefficient  of  thermal  expansion, 
a , and  a constant  coefficient  of  moisture  hygroscopic 
expansion,  3 , the  constitutive  relations  may  be  written 
for  the  lamina  in  the  material  principal  coordinates: 


el~alAT-($lra 

S11 

S12 

S13 

0 

0 

0 “ 

V 

e2-a2AT-B2m 

S12 

S22 

S23 

0 

0 

0 

°2 

e-j-a^AT-g^m 

S13 

S23 

S33 

0 

0 

0 

°3 

e23 

0 

0 

0 

S44 

0 

0 

T23 

e13 

0 

0 

0 

0 

S55 

0 

T13 

e12 

0 

0 

0 

0 

0 

S66 

T12 

k k k 

(2.1.1) 


where  the  are  components  of  the  strain  tensor.  Note 

that  the  thermal  and  moisture  expansion  effects  are  purely 
dilatational  in  the  material  coordinate  system.  The  terms 
of  the  above  flexibility  matrix  [S]  are  given  by  the  follow- 


ing  definitions:  (Please  note  the  definition  of 


’44' 


’55 


and 


S66-> 


s 

s 

s 


11 

22 

33 


1/E11 

1/E22 

!/e33 


S 


12 


~V21/'E22 


S13  'V31/E33 


S23  ~ ~V32/E33 


S44  1//2G23 
S 55  = 1/2G13 
S66  = 1/2G12 


(2.1.2) 


The  symmetry  of  the  [S]  matrix  implies  that 


v . 

l 


E.  . E.  , 
H 33 


(2.1.3) 


so  that  there  are  nine  independent  constants  for  each 
lamina.  This  corresponds  to  two  (or,  automatically,  three) 
orthogonal  planes  of  symmetry  in  the  material  coordinate 
system  of  the  lamina.  [Ref.  4) 


_ . . 


If  we  form  the  inverse  of  [S] , 


IQ] 


[S] 


-1 


(2.1.4) 


we  have 


"al~ 

a2 

°3 

= 

t23 

T13 

T 12 

— m 

k 

where 


Q 


Q11 

Q12 

Q13 

°12 

°22 

°23 

Q13 

Q23 

°33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

!11 

(1-v 

23v3; 

'22 

(1-v 

13V3. 

'33  = 

(1-v 

12v2! 

12  = 

(V21 

+ V, 

13  = 

(V31 

+ V 

'23 

(v32 

+ V 

0 

0 

0 

°44 

0 

0 


0 

0 

0 

0 

Q55 

0 


66 


e2-a2AT-62n' 

e^-a^AT-g^m 


23 

13 

12 

(2.1.5) 


(2.1.6) 
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I 

Q 

Q 

Q 


44 

55 

66 


2G 


23 


i 


2G 

2G 


13 

12 


where 


DET 


(1_V12V21  "V23V32  "V31V13  ~2v12V23V31) 


Having  defined  the  [S]  matrix  and  its  inverse, 

[Q] , in  material  coordinates,  we  now  turn  to  writing  the 
constitutive  relations  in  the  geometric  system.  The 
stress  and  strain  tensors  may  be  transformed  into  the 
geometric  coordinates  (x,y,z)  by  use  of  the  transformation 
matrix  [TJ^  to  produce  a simple  rotation  about  the  3 (z) 
axis.  (See  Figure  1) 


(2.1.7) 
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which  is  to  say  that 


[T(e)]"1  = [T  (-0  ) ] 
k K 


(2.1.11) 


Thus,  we  may  write: 


2 2 

ax  = M + N (*2 


Bx  - + N232 


2 2 
a = N a.  + M 

y i 2 


a = a0 
z 3 


3y  = n2b1  + m232 


= *3 


(2.1.12) 


axy  = MN(Qt1_a2) 


3 = MN(3,-30) 

xy  1 2 


By  combining  equations  (2.1.5),  (2.1.7),  and  (2.1.8) 


we  have  the  constitutive  relations  in  the  geometrical 
system: 


a 

X 

Qn 

°12 

®13 

0 

0 

Q16 

e -a  AT- 8 m 
XX  X 

o 

y 

°12 

^22 

°23 

0 

0 

Q26 

e -a  AT- 3 m 

y y y 

a 

z 

°13 

^23 

^33 

0 

0 

Q36 

e -a  AT- 8 m 
z z z 

°yz 

0 

0 

0 

Q44 

^45 

0 

eyz 

CTxz 

0 

0 

0 

Q45 

Q55 

0 

exz 

a 

Lxy_j 

j>61 

^62 

^63 

0 

0 

q66_ 

1. 

e -a  AT- 8 m 
xy  xy  xy 

k k k 

(2.1.13) 


where  the  entries  of  [Ql^,  are  given  by 
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[Q]k  = [T]k1  [Q]k  mk  (2.1.14) 

and  the  form  of  [Q]  is  that  of  a system  with  one  plane  of 
elastic  symmetry. 


We  may  also  invert  (2.1.13),  noting  that  an  [S] 
matrix  is  consistently  defined  by  either 


or 


[S] 


k 


(2.1.15) 


>s>k  * [<  isik  ‘Th 


(2.1.16) 


With  the  use  of  the  true  elastic  strain  quantities, 
the  [Q]  and  [S]  matrices  are  not  symmetric.  The 
non-symmetry  is  in  the  terms: 


Q 

Q 

Q 

S 

s 

s 


16 

— 

2Q61 

26 

= 

2°6  2 

36 

= 

2^63 

16 

= 

2^61 

26 

= 

2^62 

36 

= 

2^63 

If  we  had  chosen  to  use  the  "engineering"  shear  strains 
i.e.  if  the  strain  quantities  used  had  been 


Then  the  Q and  S matrices  would  have  been  symmetric. 


i 

J 


Returning  to  the  elasticity  strains  used  herein, 
we  proceed  to  invert  equation  (2.1.13)  and  introduce 

' 

our  plane-strain  assumption: 
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Defining  the  displacements. 


u = displacement  in  x 
v = displacement  in  y 
w = displacement  in  z 


direction 

direction 

direction 


we  specify,  consistent  with  the  plane-strain  assumption. 


3w 

9y 


0 


w = w (x,  z) 


(2.1.18) 


3u 

3y 


0 


u = u (x , z) 


(2.1.19) 


Thus,  we  write  the  strain-displacement  relations  for  the 
plane-strain  condition: 


e = 

X 

3u 

3x 

(2.1.20) 

e = 

y 

3v 

3y 

(2.1.21) 

e = 
z 

3w 

3 z 

1 0 

(2.1.22) 

eyz 

1 

2 

(2.1.23) 

exz 

1 

2 

,3u  3w. 

' 3z  3x 

(2.1.24) 

e 

xy 

1 

" 1 

<j^£>  -« 

(2.1.25) 

Equation  (2.1.22)  ihcludes  the  effect  of  transverse  normal 
strain. 
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2.2  Equilibrium  Equations 


Care  must  be  exercised  due  to  the  derivative  nature 
of  these  equations.  These  equations  are  written  for  a con- 
tinuum, so  they  will  be  subscripted  for  the  kfc^  lamina  to 
emphasize  that  the  derivatives  may  not  exist  at  the  inter- 
face between  laminae.  No  y variation  is  permitted  in 
d l) 

stresses,  (— = 0)  , for  the  plane  strain  condition. 
Thus,  the  equations  become: 


3 a 

x . 

Sx  h 


3t 


xz 


3z 


0 


(2.2.1) 


3t 


3t 


3x 


yz 


3 z 


= 0 


3t  . 


zx 


3a. 


»k  92  Ik 


= 0 


(2.2.2) 


(2.2.3) 


Note  that  if  the  material  is  specially  orthotropic,  i.e. 

if  the  geometric  and  principal  material  axes  coincide  so 

that  0 = 0°  or  90°  in  Figure  1,  then  equation  (2.2.2) 

is  made  trivial  by  the  proportional  relation  between  t , 

xy 

T and  the  strain  quantities  e , e which  are  assumed 
yz  ^ xy  yz 


zero. 
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2.3  Presumed  Form  for  Transverse  Shear  Stress  Txz»‘ 
Derivation  of  Normal  Stress  o z 

The  form  selected  for  the  transverse  shear  stress, 

t , is  that  of  a piecewise  smooth  quadratic  in  "z"  times 
X z 

an  unknown  function  in  "x" . Linear  superposition  is  assum 
ed  to  apply,  i.e.  the  effects  of  upper  and  lower  surface 
shears  may  be  added  in  linearly.  In  what  follows,  sub- 

4" 

script  "k"  will  refer  to  the  k lamina,  and  M is  the 
total  number  of  complete  laminae  below  the  midplane.  For 
the  laminate  shown  in  Figure  6,  M=2.  The  surface  shear 
stress  on  the  upper  and  lower  surfaces  are  t and  x , 
respectively,  and  <}>  is  an  unknown  function  of  x only. 


Due  to  Due  to  Due  to 


FIGURE  6 

Assumed  Form  for  Transverse  Shear  Stress,  X 

A 4 
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The  assumed  form  for  t^z  will  be  continuously 


differentiable  in  x and  piecewise  differentiable  in  z: 

(2.3.1) 


= +Tui*) -fj +j^Qs>l 


where 


W (':£*) 


(2.3.2) 

(2.3.3) 

(2.3.4) 


bj-5  is  a constant  which  insures  continuity  of  txz, 
Functions  f^  apd  f 2 have  been  selected  because  of  the  follov 
ing  properties  (See  Ref.  5,  p.  21)  ; 


1)  f x (— h/2 ) = 1 


f x (h/2)  = 0 


2)  f 2 (-h/2)  = 0 


f 2 (h/2)  = I" 


h/2 


h/2 


3)  / f.  (z)dz  = / f_(z)dz  = 0 

-h/2  1 -h/2 


Property  3 is  desirable  so  that  the  surface  tractions  do 
not  appear  in  the  integrated  shear  resultant- (See  equation 
2.5.31) 


By  definition,  we  must  have 


:xx\l  Hf>  - tl<x) 


(2.3.5) 


xzl  <T>  = Tu(x) 

XZ 1 2M+1 


(2.3.6) 
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We  proceed  to  solve  for  the  constants  brri  . 

55|k 

of  t through  the  laminate  gives : 

X z 

Txz|k+i (Hk+!)  - Txz|k(Hk+l) 


The  continuity 


(2.3.7) 


Using  (2.3.1)  with  (2.3.5)  we  see  that 

b , = 0 (2.3.8a) 

55  | 1 

Furthermore,  the  general  condition  (2.3.7)  will  reduce  to 


W55ik'°55ik+1)£3<Hk+l)  + b55ik  = b55lk+i 


Thus  a general  expression  for  k = 1 M is 


b55|k+1  * ^ £3«Hr+l>  lQs5lr-°5S|r+1J  <2-3'8bl 

and  there  is  reflective  symmetry  about  the  midsurface  due  to 
symmetric  material  and  geometric  properties. 


for  k = M + 2 2M  + 1 


b55lk  = b55l2H+2-k 


(2.3.9) 


Now  that  we  have  an  expression  for  x , we  turn  to 


computation  of  o . . Using  equilibrium  equation  (2.2.3), 

z Ik 

and  integrating  with  respect  to  z , we  have 


+d‘ 


(*) 


(2.3.10) 


fa fc  »■ 
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where  d,  (x)  is  a function  of  integration  to  be  determined 
z K 

and  / is  an  indefinite  integral.  Using  our  formula  for 
xz|k 

? 


ax  \ ^x*|k  d*  + dk(>) 


\T  I = -fill-  + )-dt,  /-*  +lf  + £\ 

v*l><  ix  v ♦ 2k  kv  ajr  U 21.  i.1/ 


(2.3.11) 


These  functions  of  integration,  d^  , depend  on  the  surface 
tractions  and  on  continuity  of  az: 


°*l- 


(x'Hk+l)  = ° 


*ik+1‘x'iW 


(2.3.12) 


-stm-i  ff-mn 

~ Qk)i )sr~° 


***«■  F§t+^37,-Tft4^  + 0I 


~ __ 
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Using  the  continuity  condition  (2.3.12)  for  the 
interface : 

°z (X,H2) 1 = °Z(X'U2)2 

or 


or 


This  points  toward  a cumulative  effect,  so  that,  for 
k = 2 ->  2M+1 


f irst 
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If  we  now  define 


a,  =0 


(2.3.13) 


(2.3.14) 


We  get 


(2.3.15) 


Utilizing  (2.3.15)  with  (2.3.11)  and  setting 
k = 2m+1  , z=h/2  for  the  upper  surface,  we  have 

U ' u Z*  C B 8U  8*}  &(g -*8U  K) 


(2.3.16) 


See  Section  3.1  for  a further  discussion  of  this  equation. 
Equation  (2.3.16)  shows  that  only  three  of  the  four  quanti- 
ties a , a , x , t are  independent;  the  fourth  (a  in 

J-i  U Jj  U 

this  derivation)  is  determined  from  the  other  three  by  the 
equilibrium  equation  (2.2.3)  and  the  continuity  condition 


(2.3.12) . 


With  equations  (2.3.1)  and  (2.3.8b),  and  (2.3.11) 


T 


20 

and  (2.3.15),  we  have  expressions  for  t and  ci  in  terms 

xz  z 

of  known  material  properties,  the  surface  shear  stresses, 
and  the  function  4>  (x)  . 


1 
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2.4  Calculation  of  In-Plane  Displacement,  u(x,z),  and  Local 

slope  av(x,z>|. 

8x  K 

• 3 w 

The  functions  u and  ^ will  now  be  computed  in 
terms  of  the  known  material  properties,  known  temperature 
and  moisture  functions,  the  surface  tractions,  and  the  func- 
tion <J>  . Using  the  plane-strain  constitutive  relation 
(2.1.17),  we  have 


e , = 0 = S..t  i + S . _ x | 

yz|,  44  yz|,  45  xz  | 


or 


' yz  I 


We  define 


Ak  = 


45  , 

q k Txz| 
& 44 

S 

-ill 

_ ' k 

s K 

44 


(2.4.1a) 


(2.4.1b) 


(2.4.2) 


Thus,  again  using  (2.1.17),  we  see  that 


exz|  (A+S55)k  Txz|, 

1 k 1 k 


Using  the  definition  (2.3.1)  for  x 
5S5|k°55|k  = 1 ’ we  have 


xz 


(2.4.3) 
, and  noting  that 


£x*L  = ( A * S«X  [ TL  + **  (*)] 

^ ( 1 + bQss\  ^ + + U Mk  } ^ (2.4.4) 
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Next,  we  rearrange  the  strain-displacement  equa- 
tion (2.1.24) , 

9u  | __  3 w | 

Sz'k  " ¥x  k + 2cXz|  (2.4.5) 

1 k 

and  substitute  (2.4.4)  into  (2.4.5) 

+s  (*+s5Sxyk  _ (2-4-6) 

Integrating  equation  (2.4.6)  we  obtain 

+2(i*5»iy1,d  +fkW  <2-4-7) 

Now  continue  to  reduce  the  constitutive  equation 
system  by  employing  the  other  two  plane-strain  assumptions. 
From  (2.1.17) , 

4°  _ 

(o' -a  AT-B  m),  = (S10o  + S„_a  + + S„,t  ), 

Xy  y y k 12  x 22  y 23  z 2 6 xy  k 

or,  for  the  k lamina 


r ? 


0 = 
y 


12 


22 


23 

b23  _ ^26 

°x  g-  °z  g xy 
b2  2 b22 


1 

:22 


(a  AT+B  m) 
S.  Y y 


(2.4.8) 


1 


Also,  we  have,  for  the  kth  lamina, 

o 


-ar..AT-3„..m  = S^o,,  + S^o,,  + Sc^0„  + 


xy  xy  xy 


'61  x 62  y 63  z 66  xy 


(2.4.9a) 


or,  substituting  (2.4.8)  for  , 


o - AT  + |3 


X»  ^ + ^£1  ^x 


+ s 


u 


~-b-  XT  - h±  XT  - S&r 

s„  * s * sT*-  t, 


s«  * * w~  L«y 

->21 


or 


xy 


(£„  & - sw) 


>]  + ^ + S£i 

r 

-In' 

> % + (X  - s„ 

L. 
^22  > 

(2.4.9b) 


Xz  T~ 


Cs«krs“) 


+ (SyAT+f^yhn)  _ + 


^ - S 


(2.4.10) 


Now  proceed  to  the  equation  which  will  yield  an 
expression  for  the  local  slope  • Using  the  equation 

for  the  normal  strain: 

Ez-Vr~Szm  - 513“<  + 523°y  + 533°z  * S36Txy  (2.4.11a) 


substitute  first  for  oy  , using  (2.4.8), 


a^-**,*  • v •* 


24 

f, -^AT-p,*  = S/3^ 

+ s IT  + ? T" 

J>33  Vj  ^ ->36  Lxy 
or 

£»-'*  AT- p,*,  = ( S,3-S23|i)^  + (V4|^ 

f^3S  ^3’jO^M  i^rAT+Pj,^  (2.4.11b) 


where  the  quantities  in  parentheses  are  often  referred  to 
as  the  "reduced"  flexibility  terms. 


1 
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x.  l- 

or,  introducing  a compact  notation,  for  the  k lamina, 


e | = (S'  a + S'  a + a’AT  + g’m), 

z | , 13  x 33  z z zk 


(2.4.11c) 


where,  including  a tracing  constant  "a"  to  follow  off-axis 
effects , 

r 

s 


S1 3 S13  " a 


g g (s36  s23  g ) 

23b12  + b22 


22 


(S  -S  -li) 
S66  b62  - } 

b22 


’12 


61  62  - 


22 


(2.4.12a) 


S33  S33 


- a 


(S  -S  -ii) 

q- 2 lb36  b23r  1 q 

!23  + !22_  - _g  f23} 

F + s (S63  S62  g > 

S22  ,-=■  =■  b26>  b22 


(S66_S62  - 


-) 


22 


(2.4.12b) 


a = 


a + 
z 


(S  -S  -^) 
b36  b23  - ’ 

S22 


(S  -26  _s  ) 

tS62  q S66 

22 


<*3.^*3 
22  22 


(S  -ii  - S ) 

S62  q b66; 

b22 


xy 


’23 


y s 


’22 


(2.4.13a) 


i 


J 


(2.4.13b) 


22 


S. 


(S 


_2i  _ s ) 

62  66' 


6 - — : — 6 
y c y 
22 


’22 


If  we  had  used  a plane-stress  analysis  (a  = t 

2 X Z 


tyz  = 0)  with  the  Kirchoff  linear  strain  hypothesis,  this 
would  have  lead  to  a laminate  effective  a,  "aL",  and  effec- 


tive 6,  "6  ",  and  thus  we  would  have  obtained 


e = ctL  AT  + BLm 

y y y 


= CXL__  AT  + m 
xy 


xy  xy 


e =0 
yz 


The  difference  between  the  two  methods  may  be  easily  accom- 


odated in  equations  (2.4.8)  and  following  by  substituting 


/ L 
(a  -a 


VI 


y jf , k 

for  -a 


) for  -a 


, (eL-3. 


y y 


) for  -B. 


/ L 
(a  -a 


xy  xy  | k 


xy 


, and  (B“  -3  . ) for  -3 


xy  xy|k 


xy 


For  any  symmetric 


laminate,  a^y'^xy 


0 


27 


Continuing  to  use  the  plane-strain  analysis  to  solve 


for 


3w 


3x  ' we  differentiate  (2.4.11)  with  respect  to  x , and 
use  (2.4.12)  and  (2.4.13) : 


£ r - r c;  ^ + ?'  -re*'  Q/ 

Jk  L 13  P*  Sjj  Px  + + I f 35  J 


(2.4.14) 


Now  differentiate  equation  (2.3.10)  with  respect  to  x ; 


and  also  use  the  equilibrium  equation 


in  terms  of  the  unknown  func 


tions  already  defined.  Substituting  (2.2.1)  and  (2.4.15) 


into  (2.4.14)  and  using  the  definition  (2.3.1)  for  t 
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By  use  of  the  strain-displacement  relation  (2.1.22) 

we  have 


,2 

ri  W I 

9x3z ' . 
k 


(2.4.17) 


Given  the  continuity  of  derivatives  of  "w"  within  the 
lamina,  we  may  use  (2.4.17)  with  (2.4.16)  and  integrate 
with  respect  to  z , 


IJ 

3 ? 

+ Sl3| 

+ s^lk 

1 8 GU 

z4 

_zlW2?-l 
<*/ alx* 

+s3; 

L( 

_ 

-?3 

GL 

_ z4 

+iF); 

iiT. 

:dx* 

+ u 

k2 

QssL 

z 

? 

* - 

•s'L 

J 

_Qsikd~I?) 

-cxx? 

-f- 

§4 

k ? 

3 . 

* L4,  T J 

1 m 

+ S33J 

„2 

dx 

f?  "3  AT 
J 2* 

dz 

•f 

”3  Y/i 

9X 

dz 

(2.4.18) 

where  ip.  is  a function  of  integration  to  be  found. 


The  functions  AT  and  m are  assumed  smooth 
throughout  the  laminate  and  assumed  to  be  polynomial  in  z , 
so  that  the  indefinite  integrals  with  respect  to  z may  be 


? + 


A-i^; 


(2.4.19) 


-X 


written  as 


5 If  d*  “f  ?f  * + 3C; 


(2.4.20) 


where  the  integral  evaluated  at  the  lower  limit,  z = 0 , 
will  be  zero  due  to  the  polynomial  in  z in  the  integral. 

Now  examine  the  midsurface,  z=0,k  = m+  l: 


2-'"(*.oVi  = = + X,  +X, 


If  we  redefine  the  functions  of  integration  to  include 
the  functions  and  X2  • we  simplify  the  form,  since  we 

do  not  need  to  solve  for  or  X2  • Thus,  set 


= Tk  + X,  + X. 


(2.4.21) 


for  k = 1+2M+1 


Thus  we  have 


1 , \ dw 0 

*m+i (x)  = a — 


(2.4.22) 


For  convenience,  we  define  the  following  functions 


/ -5—  AT  dz  = g (x,z) 

0 dx 


(2.4.23) 
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z 9 

/ m dz  = h (x,z)  (2.4.24) 

o 

where  g and  h are  continuous  everywhere  in  the  laminate. 

To  solve  for  the  ^(x)  in  the  other  laminae,  we  must 
make  an  assumption.  The  quantity  represents  a local 

slope  which  is  assumed  to  be  continuous  at  any  interface 
within  the  laminate.  If  the  slope  is  not  continuous, 
this  implies  delamination  and  void  formation.  To  build 
up  the  solution  for  the  > proceed  to  the  M+l,  M+2 
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+ Qss/|  (S3ll„„  S3,|„„)!^5U 

+(^KsL,-^Uu)$i]  0 
*(P'L  -P.'l  M+2 


Thus,  we  may  see  a general  pattern  for  the  cumulative 
effect  layer  by  layer.  Introduce  a tracing  constant 
"b"  to  determine  the  effect  of  transverse  normal  strain, 
ez  * use  the  form: 


y * W = cHT  -+  Ik  d#  + /<*  7*  + 7k 


dx‘ 

tS«f‘,**r«+£.5'+9/'>  (2-4‘25> 


where,  for  k=M+2->2M+l  , 
k-l 

v = U 
K 


J*  = 


5 OssL.,  S:j  ' J^r) 

H 5n5L- V bL ) ■ >L ■ * Qssl, 


(2.4.26) 


32 


r = M«-l  <•  h v 


(2.4.27) 


~h  = L (2.4.28) 

(k  DhT^H(v  »UM  <3  6^  4 1?  6/) 


^=^gjrH-§;u(f4-^)} 


(2.4.29) 


£k= v»iz  s (53x-s;l)h-] 

K r=M-H  C 


(2.4.30) 


(2.4.31) 


(2.4.32) 


vx> = ? <xA*..> 


:M*I 

+ 


(2.4.33) 


Now  examine  the  M+l,  M interface,  z=K.  . , to  set 

M+l 

up  similar  constants  for  k = M 1 (lower  laminae)  ; as 
stated  earlier,  the  slopes  must  be  equal 


Dx  (x,HM+1)  . 9x  tx'HM+l*  0 

M+l  M 


w 


I 


0 = Vmry H 

+ ( SJ  “Si  ^ / iW  + Unti  - linlL  - li  \ i-Xu 

' tolr-i  J [ 8 C,U  4K2  8 ) Jx2- 

+ ( 5 - ^ ( -^Mt|  _ Mm«i  _ Hh.i  .)■  Hm it  A Xy 

K 33  h*i  33|m/V  8 GK  4U2  8 J d%Z 

Hm*,  ( SM  Qs5|m+|~  ^ 

"f  ( QssL"  ^33  5ssln)  ( ’7m  ~ T^) 

. -( Hhf  +Qs5|(  (Kl-  s',l ) 

+ (s  L ) -s' l I ) &•  1 e£4- 

V i3bSS U j3  "551 M / 2 J 

+ c ?;i„.,-s;i„)hwi  ^ + (« 'L-  „ a j ^ «„0 

"+  ( A?|m)  ^ fX/ 

This  suggests  the  cumulative  effect  for  the  lower 
laminae.  Defining  quantities  as  in  equation  (2.4.25) 
for  k = M -*•  1 , 
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- sr  j*  7?  ^ ^ x? 


X*  ■ J"  olx'  r*  ■ IK  -J 

, + ^K0“  +s„  rl  + s„ru  +g^  +s„«  (2.4.25) 


where  we  again  use  the  tracing  constant  "b"  , to  determine 
the  effect  of  transverse  normal  strain.  We  define  the 
terms : 

L - b ^ j QssL, ' ^3  QsP  ( it  ' ) 

J*  r=M+i 

+ C^ji  Si-  SJ3“U  H.  ■‘•Qssl,  (’5|j|r.|-  S„)r  )V£ 

+(s;u.,-^y.)T 


^(k  ~\>  ) ^Tr  ( Qss]**  Qss|f-|)  J 


> =bgj(s„|,-s!;jj(^4-i-T-')}  <2-4-w 

-^k= biz  ^sJVs;3|r.l)(^-^-|+^)]  t2-4-37» 


(2.4.34) 


(2.4.35) 


sk=b|:i^5'|r-s,;L)(^-^)} 


(2.4.38) 


35 


f-Mti 


(2.4.39) 


k *1 

Cjt*^.,  U53i-?33k)^} 


(2.4.40) 


w*) 


r;M»l 


( @2  ^ (XA)  } 


(2.4.41) 


Note  that  V, 


is  the  only  constant  symme- 


The 


trie  about  the  midsurfuce,  x.e.  *J2M+2-k  ~ ^k 

are  anti- symme trie , i.e.  c = -e,_  • All  the  others 


r>k  , n.  / <5 
anti-symmetric. 


Jk  ' ,rk  ' "k  ' ~k 


'2M+2-k 

and  6”  are  neither  symmetric  nor 
k 


Thus,  we  have  implicit  in  our  definition  of 

a tracing  constant  "b"  which,  if  set  to  zero,  will  force 

all  . If  we  re-write  (2.4.18)  to  employ  the 

yk  dx 

definitions  (2.4.25)  thru  (2.4.41), 
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3Wi(X,2)K 

3* 


= av 
ax 


+[y,  k 

-Xl(r-£-^ 

¥ 

+l~^3!l» . Qss|«  ("t- ~frO~c,f|< ? 


ch* 

~ ¥¥ 
- clxz 


15,;  Q 


? +H* 


+ Qs5|, 


h-|  -+- 

3 


+ 


y.  f ] + £ 

ll  ? + *>]  f 


oll  + 


X 

/ , — 


%*(*) 


+ “>UM  + Bi|t  Wx,*) 


(2.4.42) 


where  the  function  gk  contains  both  temperature  and 
moisture  functions  of  integration  within  it. 


Now  integrate  equation  (2.4.22)  to  supply  the 
last  item  needed  to  find  u(x,z)  in  equation  (2.4.7)  . 
This  is  done  to  obtain  an  equation  for  u in  terms  of 
the  known  material  properties,  the  functions  <j> , tu 
and  a , and  the  known  functions  m and  AT  . 

Li 


' 


= ?af 

+r  s \ / _ ?g  +2^ti.\ 4. ^ 2~i 

[ ^33)k  [24  241,  201,*  iz)  'k*J  dx2 

+ [ ilaSaJ*  ?3  + /<*  *j  $* 

^ .¥+fc«k  £+!.*}$ 

♦[s^+^ja?  +i^x, 

+“ijt  ^ SM1**  * P?)^  ^ + §a^*) 


(2.4.43) 


In  an  effort  to  simplify  the  function  of  integra- 
tion, refer  to  the  integrals  below,  noting  the  continuity 
of  the  integrand; 


\ 2 = ^ el*  +JC3W 

j Vt/jiy*;di  =|o  kK2)o!*  + x4  M 


(2.4.44) 


(2.4.45) 
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where,  due  to  the  assumed  polynomial  form  in  z of  g 

K 

and  h^  , the  definite  integral  evaluated  at  the  lower 
limit  z=0  is  zero.  Now  combine  the  functions  of  integra- 
tion by  letting: 


Qoj  - ~ +X4wy) 


(2.4.46) 


where  the  C,  (x)  are  unknown  functions, 

K 


We  now  have  an  expression  for  u(x,z)k  for  a 
plane-strain  plate  element  under  the  combined  effects 
of  transverse,  normal,  axial,  temperature  and  thermal 
loadings.  Using  (2.4.7)  (2.4.43),  and  (2.4.46),  we  ob- 

tain the  expression 


u(v;K 


= C*W  - +[S»lk 


) 


+r (2  +2AQSS\  +(2A^Ss^  y,  * 


L 

-UlSslr  2*  -/i*  zl  ^ 

3k2  J 


(over) 
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-**w 


~ "^]k  J ~ PMk  )#  ^X'?J)C^  (2.4.47) 

where,  as  before,  the  function  of  integration,  Ck(x)  / 
contains  all  constants  of  integration.  Let 


— 1 ; t 

lk(X,z)  = -C"'2-|k  \ di  "P?|k  \ 


(2.4.48) 


and,  due  to  the  assumed  polynomial  form  in  z for  AT 


and  m , we  have 


Wx'0)  = 0 


Thus,  for  k = M+l  z=0  , we  have,  by  definition 


u(x'0)M+l  CM+1(X) 


CM+1(X)  = U°(X) 


(2.4.49) 


Proceed  to  the  M+l  , M+2  interface,  z-H^+2  • 
and  assume  continuity  of  u(x,z)  to  obtain  the  functions 
of  integration  C,  . 


u(x,HMJ_9)  - u (x,H  9)  _ = 0 
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| 

\ 


+ Qss),  (s»|„.f 

+ ( y m-  s 'a  KsL)  f - L,  (&)}$ 

-,WW8  'Ml)  + (^h„M»n)  ~7mjxAj) 

We  refer  to  our  definition  of  j*  p,  n»  n / <5# 

6,  £ , and  note  that  a subtraction  such  as  (MM+^-yM+2^ 
will  cancel  all  but  the  last  term  in  the  two 


series,  i.e.  the  r = (M+2)-l  term; 
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EK  "It,  {(  S'5Q»U  3^ 

+( 2 AQss|,~2  A^jS|rH)  H„, 


4- 


3k2 


+[  (2  h+  2 SjsX  (2  A +2  «», 

+ b[H„,^JH^1  (SuOss^'S^Qsslr.,)!  j 


(2 


E*=  fr^t  \ 9*L"S»  3k'rl~^  J' ) 

+ ( k 955  L ' ^ Qj5|>-+|)  (ZHni  yjy-  J 
+z[  (A+S„)r\>iSl-(b  + Sss)tu  ^+i3 


C 


I 

^ c . \ /-  Hrn. _ Hm  + Hr'.i  + {It'-l  ) 

- / (S35  y"5>jU./(  24  24U  20k'  H ' 

r*n*-i  L 


r,  , ->  — ' . / Hf.l  . Hr.i  _ Hr  *l  _ ti  Hrr  i 

t[H^,(s„|r"s»U)(x  + tv.  n*  5 


■)]) 


rsh-H  4. 

— hX.  + 1 ^ - i=  4 til-'- 1 Uk-)} 

XT  D7T^  ™l*  4U1  '6  5 


24b  £b  20b' 


.4.52a) 


f 4.52b) 


.4.53a) 


(2.4.53b) 
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G.-r 


r-Mti 


4* 


1 £ s»)r -~s» !^)(~k  r ' + t? " 1 1 tH 

f r ^ ^ -/  ft  ^ssUjC-r^  t+24")1 


Hh 


K-l 

21 

y~-n*i 


s' I \ ~— ' 
2 


•*  C K„  (Sul,-h>lr.,)HrJ} 


I ^JjL-.r5Jj|r)  (%’-!>  Hr.,)} 

=,C,  +7,(X,Hr„) 

~ £„,  6<A.i)J 


(2.4.56a) 


(2.4.56b) 


(2.4.57a) 


(2.4.57b) 


(2.4.58) 


111 


46 


face, 

laminae 

or 

0 - 


Employ  continuity  of  u at  the  M,  M+l  inter 

; = h , to  obtain  similar  constants  for  the  lower 
M+l 


0 


u (x,H 


M+l 


) " 
M+l 


u(x'hm+i)m 

M 


■*■(  (k  + S5s)H»“ 

- f)f44l 

h 

Z Wm3+,^) 

+ 

■4X; 

C + 

21.  T*~  / 

(~  Hm*I  + 

L z 

HMt|  + Hm*I  'j 

/■>  T*-  / 

,)Jn 

+[( s33)M+-s3j 

i )/-h^,-hL 

>lnJ[  2+  2+k 

u5 

+ MM4-/  •+ 

20  k2 

• W*»f  ^ 

7T  / 

wj  ax2 

^7% 

\ Wm»( 

il'vV  24  24), 

A, 

20k1 

;6  / 

_ ^ HH 

~ 11  iif" 

VJ  ax2 

+ [((2*2*^ 

+ ((>  + 2S55)M+>55L  Clh*2SSs)HKl)  ti»H 


(over) 


i 

i 


1 


i 
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* «>l-^l.d9s4C-f-ihf)  j <2 


i"  - M+-I 


-/  ■= 


S,3  Qss  f - S(J 


I — 


<?«!,)  ^ 


+ (2AQS5|,-2AQ„)m)H, 

+^<U^$5|r)$ 

+ ((^h^Ssl-(L'Sx\M-) 

k [Hr  J,1  H*  ( Qss)/’  5„  Qsj)v-|  jJ  > 


?Hr 


Ek  = ^"  ^ ^ Qssl  r-l 

K r*h-»i  l 


+ ( A Qss),'  ^Qss],..,)  *kl  ^ 


C 


^[H,(5;,|r-V.L)(f *ar  -r)]  ? 


.4.59b) 


.4.60a) 


.4.60b) 


•4.61a) 
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- $ ' I + tii-  + 

_ c ' i \ /nl 

>}l'-'){8  QU  41.* 


)('&*£)  ■ 

) + t[«.(s.'.L-5»IJ(?'-^)]I 


G.  = +tf+Pf’'3^ 

+ [ f a * Ss5),  ‘^+  j»sW]  )3 


(2.4.61b) 


(2.4.62a) 


(2.4.62b) 


(2.4.63a) 


(2.4.63b) 
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— K‘l 
A 


V 5«l.'5'iL-,)(|i:  *;r) 


__  k*l 


r = 5"  5/s'l  -s']  ^/Jk-LdZ  + H 3-3kJl/) 
b«  rw  lC  °‘r  *>  IT  U'J 

*[ (^hs)r- ( 


Hk  " ^:1,X  (■  l''S»lr-,)H']j 


K +1 


ikW=2  $-Hr( 

r-  Mtl(, 

+ (%,  (*.Hr)-Irl*,Hr))J 


(2.4.64a) 


(2.4.64b) 


(2.4.65a) 


(2.4.65b) 


(2.4.66) 
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Only  the  constants  are  symmetric  about  the 

mid-surface.  All  other  constants  - D,  , E , F , F , G 

k k k k k 

G,  - are  neither  symmetric  nor  antisymmetric. 


We  may  now  write  equation  (2.4.47)  for  u(x,z) 
with  all  functions  of  integration  known: 


k ' 


u(x,?)k  - + [_  + 


In  what  is  to  follow,  we  will,  for  convenience,  combine 
our  functions 


V)l/x.2>)  = Kt*)  ,0  A 
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2.5  Derivation  of  the  Integrated  Axial  Load,  Shear  and 
Moment  Resultants 

In  order  to  generate  sufficient  equations  to  obtain 
a solution  for  the  plate  element,  we  must  derive  expressions 
for  the  integrated  resultants  in  terms  of  the  surface  trac- 
tions, the  function  <J>  , and  the  known  material  properties 
and  known  temperature  and  moisture  distributions.  From  the 
constitutive  relation  (2.1.13),  using  our  plane-strain 
assumption , 

o 

^ 1 9"Jk  fir*, AT-ft,-)  - Q \(/y- « HT-J3  4 

* /k 

+ 1 (2-5-  la) 

or,  using  the  strain-displacement  relation  (2.1.20),  and 
omitting  the  k subscript 

“ Qii  &T-/3y*) 

^ ^T- ^ ^ ^-«xy bT~(l„ *)  (2.5.1b) 

Use  equation  (2.4.11c)  for  e and  subtract  (a  AT+8  m) 
from  both  sides 

■'***  ~ * ^}3  % 

+ (2. 4.  lid) 

The  Q term  may  be  a second  order  effect.  Introducing  a 
tracing  constant  "c"  in  order  to  determine  the  coupling 
between  axial  stress  qr  and  normal  strain  e , we  obtain 

X « 


I 


= 9«')( '“x 61  ^ +Q«|/'wr^T-^4 

'ic  *"  ^ 

* cQ.3^  [?.jL  +^»],5 

+ ft' t - ftt//-.,  (2.5.1c) 

where  az  may  be  obtained  from  (2.3.11)  and  (2.3.15), 
and  from  differentiating  (2.4.67), 


Lu—L 
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+ 


“ Q»l«  *xlr  at  - Q»L  ft L ^ 

- °^v|K  AT 

c Q>i\ 

+ £a£,+0‘J  ^ ^*t--w*lk)AT 

“ Q (&|k  wxy)k  AT  - Qisln  /^xy|k  v" 


4- 


(2.5.2a) 


Note  that  the  secondary  effect  whereby  c>x  appears  on  the 


right  hand  side  of  equation  (2.5.2)  is  cancelled  if  c = 0. 
Grouping  common  terms  in  the  above  expression  gives: 


C%X)*  I_c0 


CQ|J5u)k 


Q»L 


r40 

+ Qiilk SK  ~cQ„SuJ(  f ji ■* 'p *« jjj" 

oXU^)*U^r‘Kt  M 


(over) 


j 


p 


r 

f W 

- 

i 

i 
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We  proceed  to  form  the  integrated  stress  and 

moment  resultants 

usually  defined  in  plate  theory: 

2M+1 

Hk+1 

M = M = 

X 

I 

k=l 

/ a (x , z ) z dz  (2.5.3) 

H,  x K 

k 

2M+1 

Hk+1 

II 

2 

III 

X 

2 

II  M 

I-1 

/ o (x,z)  dz  (2.5.4) 

H.  X k 

k 

2M+1  Hk+1 
Q„  = Q = Z f 
k=l  H, 


x 


T (x,z),  dz 
xz  k 


(2.5.5) 


In  the  following  integrations,  terms  marked  (^) 
are  zero  under  the  assumption  of  material  symmetry  about 
the  mid-plane.  There  are  two  ways  this  can  occur: 


n=0 , 1 , • • • < 


h/2  2n+1 

/ p z dz  = 0 for  symmetric  "p" 
-h/2  p(z)  = p(-z) 


h/2 


2n 


/ p z dz  = 0 for  anti-symmetric  "p" 
-h/2  p(z)  = -p(-z) 


Thus 


M = 


2hn  | 


H * c*i 


% 


*Ai- 


+• 


*4 


q,l, 


(over) 


. , . | . m ^^1 


- 
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(2.5.16) 


Certain  simplications  of  these  equations  are  im- 
portant: 

1)  If  b=0 , J=0 

2)  If  AT=  AT (x) , MT  = 0 
If  m=m(x)  , Mm  = 0 

3)  If  AT=AT(z),  m=m(z),  h(x)  = 0 (2.5.17) 

Next,  we  evaluate  the  integrals  for  N: 


(over) 


L VIK  2*u  /20ll  ITS 

-$?./«>  ^?J0 

. r Co  s'  i (-£4^  +ii^Ww 
[ ) Qh  hjik  (^r  2^  ^ i6  J 

+ 1 Qnlk  fk  ~ ^ Q«/* 2 ^ J 

h[  $ 5„|,  (2t2AQ„X  (;^^)cii 

+ j Qi,/*  f2A  + 2 ?ss)x 

?„l«  010^3^  +!>  9W« 

-\QJ*#r*3i  +<0  '?  S&l*  (-QJ*** 
*5»k^-y^+V t54  )J*]g- 

'[  1 Q»U  ( Qjs)k'fr’’'Qs!Jk'0'-'<* §- 

+5U  + 

Q,,l«  *£,<<*]  3$ 

f[  -$  Q.Sj'i  +^»  h*  ^ 

-[ (qaL  +3,!“/|«-c^lK'“A+'5«“« ''‘)flTo,?  tovec’ 
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- $ (Q„  U +Q,z  (U ' c' 04  Ifi'rhl + 0.  fU, ) ■»  * 

Again  employing  simplified  notation  and  evaluating  the 
integrals  gives: 


N = 4 


ili  + 8 di-  +c  di^ 

a'/  dx3  o/x 


oU3 

^ + 0 "~T-  ' p. 

-Nt-/T  + E tr  4- U*  1 >o 


+ C + P ^ + /4  iiie 

dx 


where 


/»  = 


Z M + l 


k-i 


S„|k  *3)kL  * 6 

+ ,HA.r«0 

- Qiik  k 


(2.5.18) 


(2.5.19) 


2M*| 

H 

k=l 
+ C 


C — 


^ i-c ^js;,)k ) Q"l*Ek 

0,3  "S33U  ^ ~3  ^ssl,) 


(2.5.20) 
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•2M*I 


E = 


k = l 


l'CQ„S;3}k  l C Q'3S/ j)*  (H^"H 


(2.5.29) 


(2.5.30) 


Certain  simplications  of  these  equations  are  important. 

1)  If  c = 0 , E = 0 

2)  If  AT  = AT(x)  , m = m(x)  , then  h(x)  =0  . 
Finally,  to  determine  Q , use  the  form  (2.3.1)  for  Txz| 

k 

with  (2.3.2),  (2.3.3),  and  (2.3.4)  inequation  (2.5.5): 


or 

Q-Kf 


(2.5.31) 


where 
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+ ^SsJjc  ^fc*r  ^k)  3 (2.5 

The  simplified  form  of  equation  (2.5.31)  is  due  to  the 
convenient  forms  chosen  for  f^(z),  f ^ ( z ) ; e.g. 


2M+I 


k=l 


= 0 


«-/* 


. 32) 
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2.6  Derivation  of  Equilibrium  Equations  for  Integrated 
Resultants 

At  this  juncture,  we  will  derive  the  equilibrium 
equations  for  the  plate  element  subjected  to  surface  shear 
and  normal  stresses.  This  is  done  to  obtain  equations 
relating  the  integrated  stress  couple  and  resultants  to 
the  surface  tractions. 


> Nx  + dNx 


Mx  + dMx 


FIGURE  7 


Differential  Force  and  Moment  Components 


Assuming  a static  loading  condition. 


Fy  = 0 = £ftx  + d N y ~Tl  dx  +T'li  dxj  dy 

4-w‘  -?L+tu  = 0 


(2.6.1) 


0 =[Mx-Mx-dri,+Qxolx+Mx^ 
-V„d,&L  -Tu<Ax  j]oty 


37’  -Q*  + j(?L+tu)~0 


(2.6.3) 


These  equations  are  identical  to  those  in  [Ref.  5]  , p.  10, 

with  all  terms  relating  to  of  stresses  set  to 

zero  pursuant  to  the  plane-strain  approximation. 

Should  the  loads  be  compressive  so  that  large  de- 
flection or  buckling  results,  the  normally  second  order 
effect  of  the  vertical  component  of  the  in-plane  force 
resultant  = N becomes  important.  The  pertinent  equa- 

tion in  this  case  is: 
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X Fj  = 0 = [-  Qx  + Q,  + jqx  -v:  J*  -1  VuJx 


+ N,  -i- T dx  + N.  Is*  + 2H* 


or 


v2  o 

w 

’*^x2 

M >c 
^ J 


3x  ^x  3x 


7 
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2.7  Coefficients  for  a Specific  Hygrothermal  Case 

We  will  investigate  the  case  for  which  the  thermal 
and  moisture  effects  may  be  expressed  in  the  following  forms, 
which  are  assumed  to  be  sufficiently  general  for  all  practi- 
cal cases  of  interest. 

KtfrV  ~ + <\A  2^)  (2.7.1) 

Z r5  ■?  +/j  (2.7.2) 

In  all  the  equations  which  follow,  the  original  definition 
is  given  first,  with  the  equation  number,  and  is  immediate- 
ly followed  by  the  simplified  form  which  the  definition 
takes  if  we  use  (2.7.1)  and  (2.7.2). 

The  functions  defined  in  equations  (2.4.23)  and 
(2.4.24)  take  the  form 


i 

yo 

= J\  _(k  ? + A + g\  -As?1  + ) (2.4.23) 

' 3 2.  3 <f  ^ 

^ Is*- 0,2 

(2.4.24) 
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Now,  proceed  to  the  next  step,  equations  (2.4.27)  and 
(2.4.41).  For  k = M+  2-*-2M+l  , equation  (2.4.33) 
applies : 

k-'  i f 

-k  (2-4-33) 


For  k = M ->  1 , equation  (2.4.41)  applies: 


Now  examine  the  form  for  ^(x,  z)  , equation 
(2.4.48)  , for  k = 1 -*■  2M  + 1 


J -Z  ' f ~ 

Z^>=-°u  I j. 


(2.4.48) 


*««s"4  -A. 

-p',1 

1 3,«  ^ z s / 2 y 


In  the  solution  for  the  functions  of  integration  Ck(x) 
we  must  solve  for  £^(x) : 


(a)  for  k = M+  2->-2M+l 
k-»  c 

n ' r«  ttti  C 

(-^r  £ X,^r„)  ~ 4+1  J (2.4.58) 


{„«=},=  z:  5-H„  (3,  -J.„)  + 

* r*n*f  6 

7-  c* 7 ] +e^j  \/[  A* 4.AjHr*i) 

(.  * Iv  2lr*,J  *[  2 £ /2  20  7 

+6Pii 
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or 


(b)  for  k = M -*•  1 


k*l 


^ L ^ = l H'  (fa  I*)  -$,(*)) 

+ (^-r-'C*J  Hr)~K  (XjHr))  ] (2.4.66) 

k±I 

- ^ + ( ^ ( ?r-|  ~ %r) 


The  next  step  combines  a number  of  the  above 
functions  for  convenience: 


fyx  C*)  ? ^ £ k (t-;*)  + d k (xj  = h k (X;  2j  (2.4.68) 
or 

“*  "-Vo  + i, 

Since,  for  these  temperature  and  moisture  profiles, 

3 h 

k = 0,  we  have  no  h term  in  the  M equation  (2.5.6) 
and  no  h*  term  in  the  N equation  (2.5.18) . We  do,  however, 
have  M^,  Mm,  NT,  and  Nm  terms  shown  as  follows. 
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T 

For  M , use  equation  (2.5.15),  insert  the  AT 
form  (2.7.1),  and  note  that  the  term  (A.^  + A2*)  may  be 
pulled  through  the  summation  sign. 


z«*'  . r_ 

(a.  a,x)S  j^j,  l [q^l 

+ On  ^/)k~  C Q + P*  *xy  ]k  J * 


Similarly, 


7**-'  I f'  r _ 

1*=  ( r • r,x> r 7^5^  £ f Q.P 4 

4-  h „ rs  (C:«!) 4- r(  ti.rOj  j 


Now,  refer  to  the  N definition  (2.5.27)  and  again  insert 
the  AT  form  (2.7.1) , 
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2.8  Summary  of  Formulae  for  Stress;  Calculation  of 
S trains 

Once  the  basic  variables  have  been  obtained,  (see 
Section  3.1),  there  is  one  most  direct  order  in  which  to 
solve  for  the  stresses.  The  following  list  should  assist 
the  reader. 


a 

X 

From 

(2.5.2) 

a 

z 

From 

(2.3.11)  with  (2.3.15),  (2.3.14),  (2.3.13) 

Txz 

From 

(2.3.1)  with  (2.3.2),  (2.3.3),  (2.3.4),  and 

(2.3.8) 

Tyz 

From 

txz  , with  (2.4.1b) 

Txy 

From 

a and  a , with  (2.4.10) 

X 2 

CT 

y 

From 

a , a , r with  (2.  4. 8) 

x z xy 

Once  these  stresses  have  been  computed,  we  may 
calculate  the  strains  by  use  of  the  matrix  equation  (2.1.17) 
in  the  following  form: 


CHAPTER  3 


Studies  of  a Single  Plate  Element 


3.1  Solution  for  a Single  Plate  Element  Subject  to  Known 
Surface  Stresses 

We  examine  all  of  the  equations  and  unknown 
functions  present  at  this  step.  The  temperature  and  mois- 
ture effects  are  the  result  of  known  functions.  Assume  that 
all  surface  normal  and  shear  stresses  are  known.  There  are: 

. t m T _n\ 

Eight  known  functions:  t , t , o , a , N , N , M , m ; 

Li  U Li  U 

Six  Unknown  functions:  u°  , w°  , <J> , N,  M,  Q ; 

Equations:  3 equilibrium  equations  (2.6.1,  2.6.2,  2.6.3) 

3 definitions  of  M,N,Q  (2.5.6,  2.5.18,  2.5.31) 

6 total 


Thus,  we  seem  to-  have  exactly  as  many  unknowns  as  equations. 
However,  we  also  have  the  equation: 


o 


L 


a + 
u 


x a2M+l] 


difr 

dx 


(2.3.16) 


If  we  differentiate  the  stress  resultant  equation  (2.5.31) 


do  , 

~ , we  get 
dx 


and  solve  for 
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i 

J 


i 

j 

i 

i 

. 

* 


r 


dQ  _ . d(j> 
dx  5 dx 


and  compare  this  with  equilibrium  equation  (2.6.2), 


dQ 

dx 


o 

u 


We  can  combine  these  equations  to  get 


o-o  = kc 
L u 5 dx 


(3.1.1) 


A quick  comparison  between  (2.3.16)  and  (3.1.1)  shows  that 
we  must  have 


k5  3 °55|1  a2M+l 


(3.1.2) 


or  else  we  have  an  additional  equation  which  is  inconsistent. 
By  a check  of  examples,  including  a five-layer  laminate, 
and  by  use  of  the  symmetry  of  b^5  and  material  properties 
and  the  anti- symmetry  of  Hk  , equation  (3.1-2)  has 
been  verified.  Thus,  we  do  have  exactly  as  many  equations 
(six)  as  unknowns. 


We  will  now  proceed  to  systematically  eliminate 
variables  in  order  to  arrive  at  equations  for  the  displace- 
ment variables  u°  and  w°  . Substitute  the  integrated 
quantities  (2.5.6),  (2.5.18),  and  (2.5.31)  into  the  equi- 
librium equations  (2.6.1),  (2.6.2)  and  (2.6.3): 
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0 - A + £ ci^  + C dlfL  + V ilt*-  + c ilii 

U ^ nix1  oU*  nix  n/x2  0{  * 2 

+ n 5j7T«  + c 317-  + ~k  jzu°  _ 3 ^T-Anr 
nix1  ch*  JT  37  oTT 


0 = - IP. 


3 k* 
37 

-V 

L u 

- if. 

Ji 

1 J* 

-VZ  + t7 

0 

• + F 

+ g£* 

4-  H ill1 

ol  X4 

3*a 

3x4 

+ H 

34r, 

+ T ^ 

+ 5_k_ 

3x4 

3xz 

Jx 

y - 

" i? 

+ 1 dU 

3x2 

~ks*  *J2  (T^r“) 


(3.1.3) 


(3.1.4) 


(3.1.5) 


Solve  equation  (3.1.4)  for  . It  may  be  directly  sub- 
stituted into  (3.1.3),  and  substituted  into  (3.1.5) 
after  (3.1.5)  has  been  differentiated  once. 


£■ 1"-"' 


(3.1.4b) 


— A f 5 77  _ rl  3 P~1< 

ks  ( 5xJ  o( X ^ 


K 


3 it  -3_ra 7 +r  3l<t 
57  5x  J AIT* 


+ n j 1 ( L + r 5_1/  U 4-  b d Z ( M -f  £ 5 3 17 

5x2  5x4  5x2  5xJ 

+ y)  - 5 N T - 5 M 4 £ 37  4_ 7k  -77  +?- 

3xz  5x  o/x  5x  5x 


(3.1.6a) 


.1 
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or,  using  an  operator  notation  where  D = — 

x dx 


ilu-° 

dxa 


jL+1 

k<3  A 


-[ 

+4  r 4NT+  cUC  _ Jj*  1 

d X dx  J 


(3.1.6b) 


Now  proceed  to  equation  (3.1.5),  differentiating  once  and 

then  substituting  (3.1.4b)  for  ^ 

dx 


= -p„  4V  + JVJ 

dx4-  ^ [d?  J7  J 

+ ii.  d Vu  _ olZUl~j  + H d 5r,  4 1 ddjl 

ks  Ldx2  dxz  J dx-5  (dx1 


+ M -f 

dx5 


T dlru  +dM,  _ dVf 

dx3  dx2  dx2 


or,  again  using  the  operator  notation  D = — 

x dx  ' 
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+ -i- 


ci'T:  _<j 


t),,  L A: 


<jU2  d-  x2 


(3.1.7b) 


Thus,  with  the  right  hand  sides  of  both  (3.1.6b) 
and  (3.1.7b)  known,  the  solutions  in  terms  of  displacements 
are  given  by 


x x 


u° (x)  = Ag  + A^x  + / / right  hand  side  of  (3.1.6b)  (3.1.8) 

P - x x x x 

w°  (x)  = Bq  f B^x  + B2 x^  + + f f f f right  hand  side 


of  (3.1.7b) 


(3.1.9) 


x 

where  / indicate  an  indefinite  integral. 


The  evaluation  of  these  constants  in  terms  of  real  boundary 
condition  will  not  only  involve  w°  and  u°  , but  also  Q 
M,  and  N.  It  is  necessary,  therefore,  to  derive  expres- 
sions for  M,  Q,  and  N in  terms  of  u°  and  w° . 


rr 

in> 
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(3.1.4b)  for 


Return  to  equation  (2.5.6)  and  substituting  in 
d4> 


dx 


+ +I 
+ H 7s“+Jzr 


(3.1.10a) 


or;  with  the  operator  notation  Dx  = JJx 


h = 


+■ 


n pj  + /£ 


rL 


+|7/op*  +li-nr-ti~  ] 


(3.1.10b) 


A similar  substitution  may  be  performed  with  equations 
(2.5.18)  and  (3.1.4b) 


fill +1 

olx*  dxz 


(vL-\ru) 


n = A 

Ks 

+ C + P + C -t-p  dXu 

c[xJ  %x  ^ dx 

+ E C-J-  + A s[*  - Nt-  r+  Vi”f0  + E o: 


(3.1.11a) 


or, 


jl 
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N given  by  (3.1.11b) 

M given  by  (3.1.10b) 

Q given  by  (3.1.12b),  or,  more 

simply,  by  integrating  (2.6.2) 


and,  rearranging  (2.5.31), 


(2.5.31b) 
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3.2  Parametric  Studies  of  a Single  Plate  Element  Subjected 
to  Known  Surface  Tractions  and  Known  Temperature  and 
Moisture  Profiles 


In  the  interest  of  demonstrating  the  practicality 
and  consistency  of  this  method,  several  typical  examples 
have  been  selected.  These  examples  include  a plate  sub- 
jected to  simple  normal  tension,  simple  shear  stress,  and 
various  temperature  and  moisture  profiles.  Where  a specific 
example  is  needed,  material  properties  for  a laminate  of 
T300/5208  graphite  epoxy  are  used.  (See  Fig.  8)  Figures 
9 and  10  will  help  visualization  of  the  boundary  conditions. 

The  heights  and  laminae  are  numbered  as  shown  in  Fig.  5; 

-3  . 

h.g  = .1397  x 10  m is  the  thickness  of  each  lamina. 


FIGURE  8 


H6  - . 4191 x 10'3  m 
H5  - . 2794  x 10' 3 m 
H4  * . 1397  x 10'3  m 

H3  - 1397  x Kf 3 m 

H2  - 2794  x 10"3  m 

Hj  • 4191  x 10" 3 m 


Example  Symmetric  Laminate 


FIGURE  10 


Cantilever  Plate 

Example  1:  Simply  supported  plate,  surface  shear  xL  = tu 

t . 

This  example  is  chosen  to  verify  that  the  results  of  this 
method  are  consistent  with  physical  intuition.  The  dis- 
placement u = u(z)  and  the  shear  stress  t will  be  cal- 
culated . 
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= 0 


The  boundary  conditions  are: 


u°(0)  = u°(L)  = 0 

w°  ( 0)  = w°(L)  - M (0)  = M (L) 


By  3.1.9  and  3.1.10b, 

w° (x)  = 0 

By  3.1.8,  with  t = tt  = t 

U 1_J 

u°(x)  = o 

By  2.4.67,  with  AT  = m = 0 

u / ^ = f ('I 'f  + 


+ ^2  + 2 L Qss\  2^)  + ft  ^ ^Si)k U 

- **  + *]  f 


where,  by  2.5.31b  and  3.1.12b, 
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T 


U IT 

(10"14m/Pa) 


n-,2 


FIGURE  II 

Simply  Supported  Plate,  All  0°  Layers 
Plate,  Loading:  Surface  Shear  Stress 

This  practical  example  shows  that  the  linear  terms  in  z are 
by  far  the  leading  terms;  deviation  from  linearity  in  u is 


<1%  . With  txz  = Q44*/  + Q55exZ  = Q55  \ + J$>  ' Txz 

is  virtually  constant,  as  expected. 

Example  2;  Simply  supported  plate,  normal  tension  °u  = aL  = 
<T  . This  example  is  chosen  to  check  the  resuxts  of  this 
method  with  physical  intuition.  We  will  compute  w=w(z)  and 


also  N 


u° (0)  = u°(L)  = 0 

w° (0)  = w°(L)  = M (0)  = M (L)  = 0 


AT  = m = 0 


* • i 


I 
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By  3.1.6b  and  3.1.8, 


u (x)  = 0 


By  3.1.7b  and  3.1.9, 


w (x)  = 0 


By  2.4.67, 


u(x)  = 0 


By  3.1.12b  and  2.5.31b, 


Q = <t>  = 0 


Therefore,  by  2.4.42, 


3w 

3x 


= 0 , 


or 


w = w (z) 


,0 


For  the  case  of  all  0 layers  of  identical  material. 


w (z)  = w + / ez  dz 

o 


= / e dz 

z 
o 


By  2.3.11  and  2.3.15 


o2  = a 
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By  2.5.2b, 


°x  C Q13S33°L  C Q13S33° 


h/2 

Thus,  N = / a dz  = h o > 0 , which  is 
-h/2  X X 

as  expected;  the  plate  must  be  in  tension  to  stay  in  place. 


Finally,  by  2.4.11c, 


e 


z 


e 


z 


S13(C  Q13S330)  + S33a 
S^3(l+c  Q13Si3)a 


Therefore, 

w(z)  = e z z = (j:oz  ; (j:  > 0 

as  expected.  For  T300/5208,  this  relationship,  with 
tracing  constant  c = 1 , becomes 

w ( z ) = (8.95  x 10_11/Pa)  a z 

and  for  c = 0 

w(z)  = (8.83  x 10_11/Pa)  a z 


Example  3:  Residual  Thermal  Stresses 

For  many  advanced  composites,  the  stress  free  tem- 
perature is  an  elevated  cure  temperature,  and  a large  nega- 
tive AT  loading  is  present  at  room  temperature.  For 


residual  thermal  stress  calculation,  we  may  visualize  the 
plate  as  a cantilever  (see  Fig.  10).  Stresses  are  generat- 
ed within  the  plate  due  to  the  mismatch  in  thermal  expansion 
coefficients  between  layers. 

As  boundary  conditions,  use: 

w0(°)  = gjr  (0)  = u°(0)  = 0 

M (L)  = N (L)  = 0 

AT  = constant  , m = 0 

mT  = M™  = N™  = 0 , NT  * 0 


By  3.1.3  and  3.1.11a 

T 

u (x)  = — x 
A 

By  3.1.9,  3.1.12b,  and  3.1.10a 
w*"*  (x)  H 0 

By  3.1.12b  and  2.5.31b, 

<P  = Q = 0 


By  2.5.2b, 


(Qm^xIk  + Qjz'*y)k 
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Example  4:  Structural  plate  element  subjected  to 
temperature  increase  and  moisture  in- 
filtration. 


For  this  case,  the  plate  is  constrained  by  the  sup- 
ports and  then  uniformly  heated  and  non-uniformly  infiltrat- 
ed by  moisture.  This  example  is  chosen  to  illustrate  an 
application  of  current  interest:  hygrothermal  effects. 


Boundary  conditions: 

u°  (0)  = u°(L)  = 0 
w°(0)  = w°(L)  = 0 
M ( 0)  = M(L)  = 0 


AT  = constant  = 8 0°C,  m = m(z) 
MT=0  NT  * 0 M1"  j-  0 Nm^0 


Proceeding  as  in  example  3,  we  find 


u (x) 

0 , , 
w (x) 


= 0 


= ( 


M™!, 

2D11 


) x - ( 


M™ 


2D 


) x' 


11 


4>  = Q = 0 


By  2.5.2b, 

- f AT 

~(Q"  PxL  w] 


« 4 


FIGURE  14 


Structural  Moisture  and  Thermal  Stresses 

If  a plate  is  cooled  from  curing  temperature, 
placed  in  a structure,  and  subjected  to  moisture  and  ele- 
vated temperature,  the  results  can  be  analyzed  by  super- 
imposing examples  three  and  four. 

For  cooldown  from  a cure  temperature  of  350°F, 
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the  residual  curing  stresses  are  the  same  order  as  the 
structural  stress  due  to  moisture  infiltration.  Any 
analysis  which  ignores  the  moisture  effect  is  thus  poten^ 
tially  hazardous. 


i 


CHAPTER  4 


Derivation  of  the  Governing  Equations 
for  a Bonded  Joint 


4.1  Development  of  Alternate  Expressions  Employing  o^ 

In  preparation  for  a discussion  of  bonded  elements, 
we  need  to  develop  expressions  written  in  terms  cu»  t , 
x rather  than  oT  , x , xT  . This  may  be  done  by  employing 

L Li  U Li 

(3.1.1) 


. d <J> 

Or  - 0 = kc 

L u 5 dx 


(3.1.1a) 


or 


. , d* 

Or  = 0 + -rz- 

L u 5 dx 


(3.1.1b) 

Examine  equation  (2.4.67) , taking  first  the  pertinent  term 
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Similarly,  we  may  take  equation  (2.4.42)  and  alter 
its  form  by  using  (3.1.1)  as  follows: 


”3w(x,J)K  = cUi° 
oU 


l >i3l*  ( ¥ U H1  sj  7*J 
/ ^*/*(  S U 8/ *>J  *37» 


(over) 
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+ Wi  I + ^s  5,31,  **!</.,  5 
+ Vllj]  (4-1'2> 


Next,  examine  equation  (2.5.18): 


or,  letting  A*  = A + E 

B*  = B + E k5 


(4.1.3) 

(4.1.4) 


*& 

A *’ 


\*4£+r££r  4 1)^1 

d X 


N=  /T  VS  4 B’f^C  +/> 


4C  4l£»  + D 4 f ^ -i* 


d.*’ 


olx 


elX 

d* x 


+ /)  - NT-  hi”  + E 

£> IX 


(4.1.5) 


Equation  (2.5.31)  remains  unchanged,  containing  no  explicit 

a . Equation  (2.5.6)  also  remains  unchanged. 

Li 


M- 


3 


4.2  Development  of  Interface  Equations  for  Bonded  Joints 

The  fundamental  assumption  on  interaction  between 
bonded  elements  is  that  the  adhesive  is  a thin,  homogeneous 
isotropic  material  in  a state  of  simple  shear  and  normal 
stress.  Examine  two  elements,  "i"  and  "j"  , shown  in 
Figs.  15  and  16  with  positive  stress  and  displacements. 


Eploded  View  of  Bonded  Joint 
Including  Adhesive  Stresses 
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u. 

1 


(4.2.2) 


wl^ere  the  constants  Cg,  C N must  be  determined  experimentally 
in  thin  film  configuration  (not  in  bulk  form)  and  are  pro- 
portional to  the  film  thickness. 

We  now  proceed  to  develop  equations  (4.2.1)  and 
(4.2.2)  in  terms  of  the  unknown  functions.  First,  differ- 
entiate (4.2.1)  with  respect  to  x , since  we  have  ex- 

• 9 w 

pressions  for  ^ , (2.4.42)  and  (4.1.2). 


dx 


.\  r i_ 

' C N 


iv, 

cL* 


(4.2.3) 


Evaluating  (2.4.42)  for  element  "i"  : z 


= d w;/x  -ki\ 

I k«i 


oTx 


4 


f 5,,|,  . 3^./,) 

L ^ li 


dw;  l*/\)  — 


ct< 


+ 4s 


4 +*,]?.  i+s ^sjtK! 
Jip.‘  s' I +•,]■&  S-jJW  ua'ldb..- 

L 192  it/,  7‘J  o/x1  1 132  «/,  ^'J  dS 

+[  5,'j  ?«],  +/<,] 

+ I'ttf&H*  J.  1 ff 


+ (r  s^'+f)af  *>« 
+“»)jcvV;+  /U,I<V  kj 


(4.2.4) 


Similarly,  we  may  use  expression  (4.1.2)  for  element  "j 


with  k = 2M  +1,  z = +h./2  . Note  that  in  the  term  for 
2 1 

— , we  will  use  symmetry  and  substitute  Q^ci  for 

dx  bbl2M+l 

| . This  is  done  to  avoid  confusion,  so  that  all 

quantities  subscripted  "1"  denote  element  "i",  and  all 

quantities  subscripted  " 2m+1"  denote  element  "j"  . 


I k-2Mtl 


(*,  r) 


I *»hltn  2 J oT^ 


+-  c><. 


2«tl 


(4.2.5) 


Substituting  equations  (4.2.4)  and  (4.2.5)  into  equation 
(4.2.3)  gives: 


CN  oU  d X 


4[  Su  <?«),  f Mj  h 4[" ^ 9s5)i  *f']  J? 

+ 4 5^-tO-AJ,  Wry 


. f *7  J 
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or,  grouping  terms  for  convenience, 


where 

K --  -f  Kl ,*s,  -rKu-L., 

7j  [32  *”l  "'7»  ' 737  smL,“ 

~ Qssl  + f1! 

^5  r ^ 9«)i  *// 


(4.2.6) 


(4.2.7) 


(4.2.8) 


(4.2.9) 


(4.2.10) 


(4.2.11) 


' 
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We  will  now  perform  similar  operations  on  the  equa- 
tions for  u to  obtain  the  equation  (4.2.2)  in  terms  of 
our  unknown  functions.  Using  equation  (2.4.67) 


k-  I 


Next,  use  equation  (4.1.1)  with  k = 2M  + 1 , z - h^/2  . 

We  use  symmetry,  Qcci  = Qcc  i • t^ie  term  multiplying 

2 bbll  d:jI2M+1 

^ ^i  to  avoid  confusion. 


1 

i 
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2 


(over) 
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Ki 

4 


i^  + Sss), 


+k 


'IMtl 


(4.2.31) 


= ’420_S33^m(  F^+2  7*»'i 


(4.2.32) 


~?\i  (*)  r ^ i (Xj'^A)  ~ i &A/0 


(4.2.33) 
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4.3  Discussion  of  Bonded  Joint  Equations 

We  now  summarize  the  unknown  functions  in  a bonded 
joint.  For  elements  "i"  and  "j"  , (Fig.  15)  we  have  14 
unknown  functions . 

u°,  u°.  w°.  w N.,  N.,  M.,  M.,  Q.,  Qj.  T0,  o0 
and  six  known  functions 


„,T  „m  WT  ,rn 
N , N , M , JVT  , T . , T • 

U / 1 Lt  f J 


The  equations  for  this  system  are: 

3 equilibrium,  element  i ’'j 
3 equilibrium,  element  j J 
3 definitions  of  , N^, 

3 definitions  of  , N ^ , Q.. 


(2.6.1,  2.6.2,  2.6.3) 


(2.5.6,  2.5.18 
2.5.31) 


2 interface  equations  (4.2.6,  4.2.21) 
14  total 


W.J.  Renton  [Ref.  8]  has  found  that  the  reduction 
may  be  accomplished  by  obtaining  two  linear,  constant  co- 
efficient, ordinary  differential  equations  in  the  two  un- 
knowns Tq  and  . In  fact,  either  or  may  be 

used,  since,  by  the  EFy  equilibrium  equations  for  the 
adherend  elements  i and  j , (See  Figs.  15  and  16) 
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and 


0 


°0  + 


o 

u,i 


°T  • + 

L , 3 


yielding 


dQi 

dx  dx  + °L,j  °u,i 


If  the  adherends  "i"  and  "j"  are  identical,  the  equa- 
tions uncouple.  For  dissimilar  adherends,  using  operator 
multiplication  will  reduce  the  system  to  one  equation  in 
one  unknown,  an  eighth  order  ordinary  differential  equation 
in  t Q . 

Unknown  constants  in  the  system  can  be  evaluated  as 
followed:  3 unique  constants  for  each  element  (match  either 

u°,  w^ , and  , or,  preferably,  M,  N,  and  Q as  boundary 

conditions  at  each  end  of  the  element) , and  (x=0)  = 

(x=L)  = 0 where  L = length  (x  distance)  of  the  interface. 
This  last  provision  is  necessary  to  have  be  zero  on 

each  end  of  the  bond  where  there  is  a free  surface.  The 
main  difficulty  in  the  solution  is  algebraic,  since  we  must 
derive  equations  involving  and  each  of  the  other  quan- 

tities of  interest  so  that  we  may  use  the  boundary  condi- 
tions on  the  elements  i and  j to  solve  for  the  unknown 
coefficients  in  . 
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CHAPTER  5 


Conclusions 

An  analytical  model  has  been  developed  for  the 

plane-strain  analysis  of  symmetric  anisotropic  laminated 

plates  subjected  to  a full  range  of  thermal,  moisture, 

and  mechanical  loadings  including  surface  tractions.  This 

analysis  includes  transverse  shear  strain  e and  trans- 

xz 

verse  normal  strain  ez  , and  can  be  used  for  any  laminate 
which  through  engineering  judgement  can  be  assumed  to  be 
in  plane-strain  in  the  plane  of  loading.  The  methods  of 
analysis  have  been  demonstrated  to  be  accurate  for  a singl 
laminated  plate.  The  example  problems  (section  3.2)  show 
the  method  to  give  results  which  are  consistent  with  physi 
cal  intuition.  In  addition,  problems  of  current  interest, 
including  residual  curing  stresses  and  hygrothermal  load- 
ings (combined  elevated  temperature  and  moisture)  have 
been  examined. 

The  edge  loads  on  the  y = constant  edges  arise 
as  a natural  consequence  of  the  plane-strain  assumption, 
and  the  displacement  "v"  in  the  y direction  is  zero. 
Should  the  displacements  or  loadings  in  the  y direction 


I 


become  important,  then  the  plane-strain  assumption  herein 
will  no  longer  be  valid. 


The  results  of  calculations  involving  the  tracing 
constants  a,  b,  and  c show  that  their  effects  are  im- 
portant. The  difference  in  material  constants  is  greatest 
in  (the  stretching-shearing  compliance  term  corrected 

for  off-axis  effects,  equation  2.4.12a);  there  is  a factor 
of  two  difference  in  depending  on  whether  a = 0 or 

a = 1 . For  constants  a=b=c=l  or  a=b=c=0  , the  constant 

d3x 

"C"  (multiplying  L in  equation  2.5.18  for  N ) varies 

dx3 

by  a factor  of  five.  There  are  many  other  constants, 
however,  which  do  not  vary  appreciably  with  the  values 
of  a,  b,  and  c . Thus,  whether  a,  b,  and  c are  one  or 
zero,  the  stress  and  displacement  calculations  only  vary 
by  about  3 to  4%.  To  take  advantage  of  algebraic  complete- 
ness, constants  a,  b,  and  c should  always  be  set  equal 
to  unity;  the  computer  program  in  Appendix  B has  these 
constants  set  equal  to  one  as  default  values. 


There  are  many  potential  uses  of  this  plate  element 
as  an  "analytic  finite  element".  While  the  usage  of  ag- 
gregations of  these  plate  elements  is  only  suggested  in 
Chapter  4,  the  single  lap  joint  (see  Fig.  3)  as  well  as 
the  double  lap  or  step-lap  joint  may  be  analyzed  by  the 
methods  of  analysis  developed  herein.  The  only  difficulty 


i tin  nnnn^i— 
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with  such  solutions  is  the  length  of  the  algebraic  mani- 
pulations needed;  there  is  no  difficulty  with  the  concept 
of  using  the  plate  as  a finite  element  for  the  constructions 


discussed  above. 
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APPENDIX 


Comparison  of  this  analysis  with  previous 
analysis  by  W.J.  Renton  and  J.R.  Vinson  as  presented 
in  AFOSR  TR  75-0125,  Aug.  1974 

I.  Use  of  Tracing  Constants  a,  b,  c 

Several  terms  were  omitted  in  Renton's  analysis 
[Ref.  8]  but  have  been  included  in  this  report.  The  work 
in  this  report  is  algebraically  complete,  and  the  tracing 
terms  should  all  be  set  equal  to  1 to  take  advantage  of  this 
fact.  All  tracing  constants  would  be  set  equal  to  zero  to 
give  the  Renton  results.  This  difference  will  affect  most 
calculations  <5%  . See  equations  (2.4.12),  (2.4.26),  and 
(2.5.1c)  for  these  constants. 

II.  Plane-Strain  Assumption 

Renton  assumed  that 
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There  is  some  difficulty  with  this  form  of  assumption , 
however.  For  example,  with  a constant  AT,  and  dissimilar 
layers  i and  j , we  have  t (ay)  * If»  as  stated 

by  Renton 


*ey*k  ( 9y  Jk  (cV 


AT 


k = i,  j 


then  the  layers  will  have,  as  y varies, 

vt  * vj 

at  the  interface.  The  layers  will  thus  slide  freely  over 
each  other  and  delaminate. 


The  assumptions  in  analysis  herein  are  those  of 
a classical  plane-strain  condition,  including  moisture 
effects  in  the  constitutive  equations: 
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III.  Inclusion  of  Moisture  Loadings 


Renton's  analysis  makes  allowance  for  thermal 
strains  only,  with  constant  thermal  expansion  coefficients. 
This  report  includes  thermal  strains  with  constant  thermal 
expansion  coefficients  and  also  moisture  strains  with  con- 
stant moisture  expansion  coefficients. 


